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We present a theory for the coherent magnetization dynamics induced by a focused ultrafast 
laser beam in magnetic films, taking into account nonthermal (inverse Faraday effect) and thermal 
(heating) actuation. The dynamic conversion between spin waves and phonons is induced by the 
magnetoelastic coupling that allows efficient propagation of angufar momentum. The anisotropy of 
the magnetoefastic coupfing renders characteristic angle dependences of the magnetization propa¬ 
gation that are strikingly different for thermaf and nonthermaf actuation. 
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Introduction — Since the discovery of laser-induced ul¬ 
trafast spin dynamics in Nickel by Beaurepaire et al. [1], 
the spin manipulation in ferromagnetic system by fem¬ 
tosecond laser pulses has attracted much attention since 
combining the intellectual challenge of new physics with 
the application potential of ultrafast magnetization re¬ 
versal [2]. Intense light can cause many effects in mag¬ 
nets, such as the coherent inverse Faraday effect (IFE) 
as well as the excitation of the coupled electron, magnon 
and phonon subsystems on various time scales. The asso¬ 
ciated modulation of the magnetic anisotropy and mag¬ 
netization modulus allows coherent control of the mag¬ 
netic order [3, 4]. The transient magnetic field generated 
by the IFE allows non-thermal ultrafast magnetization 
control [5, 6] that may be distinguished from heating- 
induced effects by switching the light polarization. Nev¬ 
ertheless, heating is essential for light-induced magneti¬ 
zation reversal [7]. Toggle switching of the magnetization 
by heat alone has also been reported [8]. Understanding 
and controlling the relative magnitude of thermal and 
non-thermal excitation is therefore an important but un¬ 
solved issue. 

The optical ultrafast pump-probe technique as shown 
in Fig. 1 is an established powerful method to study 
matter. Here we will show that the symmetry of the 
spatiotemporal magnetization distribution excited by a 
focused laser beam reveals the relative contributions of 
thermal and non-thermal excitations. This phenomenon 
originates from the magnetoelastic coupling (MFC) [9- 
11], i.e., the coupling between spin waves (magnons) and 
acoustic lattice waves (phonons). 

In the regions of phase space in which the magnon 
and phonon dispersion come close, the MFC hybridizes 
magnons and phonons into coherently mixed quasipar¬ 
ticles [“magnon-polarons”(MBs)]. This implies that 
phonons in magnets can be converted into magnetiza¬ 
tion and become detectable via magneto-optical [12-14] 
or electrical [15, 16] techniques. In this Letter, we present 
a study of the spatial magnetization dynamics in mag¬ 
netic thin films after focused-laser excitation [17]. We 
consider here magnetic insulators that are not affected 



FIG. 1: Pump-probe study of the dynamics of magnetic films 
by pulsed lasers: An external magnetic field aligns the equi¬ 
librium magnetization along the x direction. The pump laser 
hits the sample at the origin in time and space. The tem¬ 
poral distribution of the out-of-plane (z) component of the 
magnetization is detected by the Kerr rotation of the linearly 
polarized probe pulse. 

by conduction or photo-excited free carriers. We find 
that laser-generated phonons efficiently excite magneti¬ 
zation when the diameter of the laser spot is comparable 
with the wavelength of the MBs. The spatial dynamics 
of such phonon-induced magnetization shows a different 
angular symmetry from that of the magnetization gener¬ 
ated directly by laser via the IFE. Dedicated experiments 
suggested here can therefore help understanding the fun¬ 
damental nature of light-matter interaction in magnets. 

Formalism — The basic theory for MBs in special sym¬ 
metry directions by Kittel [9] and Akhiezer et a7[10], 
extended by Schlomann [11] to arbitrary propagating di¬ 
rections, was developed more than half a century ago. 
The energy density of the minimal model reads [9] 

TL = Hex + Hz + He\ + H mec + Hdip. (1) 

We adopt a cubic unit cell and consider the thin film 
limit in which the magnetization is spatially constant 
over the film thickness. This assumption holds for films 
up to 100/rm [17], which for wide-gap insulators is still 
less than the penetration depth of the light and al¬ 
lows us to use a two-dimensional model. With ex¬ 
ternal magnetic field H and equilibrium magnetization 
vector Moljx (|Mo| = Mq saturation magnetization) 
Hex = A[(Vmy)^ -I- (Vto^)^] and Hz = -roHMo + 
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{fioHMQ/2){my+ml) represent the (linearized) exchange 
and Zeeman energies, respectively, where my and are 
the transverse magnetization components of m = M/Mq. 
'Hei is the lattice energy with both kinetic and elas¬ 
tic contributions Hei = (l/2)pR-R+ (l/2)A(^j -I- 

^ij with strain tensor Stj = [diRj + djRi)l2 and 
R representing the lattice displacement with respect to 
equilibrium, p the mass density, and A and p, elastic con¬ 
stants. The MEC in Eq. (1) reads "Hmec = j2t,je{y,z}i^ + 
adij)Sijmimj -I- SixTUi, where a and b are 

magnetoelastic coupling coefficients. By adopting the 
short-wave length limit of the magnetostatic dipolar in¬ 
teraction "Hdip ~ (pqMq /2)mySm^ 9, we disregard the 
Damon-Eshbach surface modes [18] and simplify the dis¬ 
persion of the volume modes, which is allowed for small 
laser spot sizes with response being dominated by high- 
momentum wave vectors [19]. 

By introducing the forces and torques F acting on 
the displacement vector $ = Ri, Rt, Rz)"^, one 

can write out the linearized equations of motion as 
shown in the Supplemental Material [20]. Here, the 
lattice displacement is rewritten in the form of longi¬ 
tudinal (Ri), in-plane transverse (Rt) and out-of-plane 
transverse (Rz) modes. Strictly speaking, the damp¬ 
ing of phonons and magnon are not necessarily inde¬ 
pendent, since magnetization is affected by phonon at¬ 
tenuation via the MEC [27]. We treat Gilbert damp¬ 
ing constant a and phonon relaxation time Tp as inde¬ 
pendent parameters since Gilbert damping can also be 
caused by magnetic disorder, surface roughness or de¬ 
fects [28]. We define the anisotropic spin wave frequency 

flo = IPe \JH{H + Mq sin^ 9) and the MEC frequency 

parameter A{k) = \//{AMopQo) with 9 being the 
angle between magnetic field and in-plane wave vector k. 

The spatiotemporal dynamics of $(r, t) reads 

$i(r,t) = J dr'dt'Gij{r-r',t-t')Fj{r\t'), (2) 

where Gij are the components of the Green function 
matrix (propagator) associated with the magnetoelastic 
equations of motion specified in the Supplement [20]. A 
femtosecond laser pulse generates forces via the inverse 
Earaday effect [2, 17, 29] and heating [14, 30-32] that 
are instantaneous on the scale of the lattice and magne¬ 
tization dynamics. The relative importance of these two 
mechanisms depends on the material and light and is still 
a matter of controversy. Here we find that spot excitation 
of thin magnetic films is an appropriate method to sep¬ 
arate the two, since they lead to conspicuous differences 
in the time and position dependent response. 

We consider circularly polarized light along z that by 
the IFE generates an effective magnetic field along the 
same direction. For a femtosecond Gaussian pulse with 
spot size W, the generated magnetic field has a spatial 
distribution HiFE(r,i) = ^^^iFE/(i) exp(—r^/lE^) where 


temporal shape f{t) « TiS{t) with pulse duration rp, the 
amplitude i7iFE = PlinO' is proportional to laser intensity 
(/in) and IFE coefficient (/?), respectively. The latter is 
related to the Verdet constant (V) as /3 = EAo/(27rco), 
where Aq and cq are wavelength and velocity of the 
light [33]. a = 1(—1) for left(right)-handed polarization. 
The torque F„j^(r,t) = -fTipoHwE.5{t) /W"^). 

On the other hand, the light pulse generates a sud¬ 
den increase of the local lattice temperature ST{r,t) = 
[r/inT//(pC'„)]0(t) exp(—r^/IT^), where T is the light- 
absorption coefficient. By choosing the Heaviside step 
function 0(t) we assume that the lattice locally equili¬ 
brates much faster than the response time of the coherent 
magnetization (a few picoseconds [34]), while the subse- 
qent cooling of the lattice by diffusion is slow. The re¬ 
sulting in-plane thermoelastic stress = rj{2>X+2p,)drT 
generates longitudinal (pressure) waves [31, 35], where p 
is the thermoelastic expansion coefficient. The local ther¬ 
mal expansion also generates a “bulge” shear stress [36] 
at a free surface, i.e., an out-of-plane displacement Rz- 
Fji^ = CryV^T, where C is a parameter proportional to 
the film thickness and controlled by the substrate and an 
eventual cap layer, leads to displacement proportional to 
local temperature gradient (see numerical results below). 
In the Supplemental Material [20], we specify the mate¬ 
rial parameters for yttrium iron garnet (YIG) adopted in 
our calculations. 

One-dimensional dynamics — We start with a spin 
wave propagating along the external magnetic field, i.e., 
9 = 0, which by symmetry couples only with the trans¬ 
verse phonons. The IFE generates the torque F^y {x, t) = 
moS{t) exp(—x^/IT^) with toq = yr/yxoi/iFE- This can be 
realized by a line-shaped excitation spot [17]. 

The calculated magnetization profiles at Hot = 50 and 
75 without and with MEG are plotted in Fig. 2(a) and (b) 
separately for a = 10“^. Without MEG the magnetiza¬ 
tion is localized at the exposure spot and broadens only 
very weakly with time while the MEG strongly enhances 
the broadening of primary magnetization packet, with a 
wavefront propagating with the sound velocity Cj. This 
phenomenon illustrates that the lattice plays an essential 
role for spin transport in magnetic films. 

Fig. 2(c) illustrates the sound-assisted propagation for 
Gilbert damping a. = 0.1. Instead of the expanding wave 
front in Fig. 2(b), we now find two packets escaping the 
excitation region into opposite directions. The packets 
have a much longer lifetime than the coherently gener¬ 
ated IFE magnetization, hence dominate at long time 
scales. This behavior is recovered by the asymptotic 
expression obtained when a 3> A(fcc) at the magnon- 
phonon dispersion crossing wave vector {kc = ^/ct), 

mz{x,i) ~ moe~°*~^ sin (t) -|- (2a)“^TOoe“*^^*’A^(fcc) 

2ac^ [Ai{cti - x)Ai{ctix)] , c* < 1, 

[A2{cti-x) -I- A2{cti+x)] , c* > 1, 
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FIG. 2: (Color online) One-dimensional model for the dynam¬ 
ics of the out-of-plane magnetization component induced 
by the inverse Faraday effect: (a) At times Qot = 50 and 75 for 
a Gaussian laser intensity spot in the absence of MEG with 
Gilbert damping a = lO”"*. The dashed envelopes are the 
modulus (|in — ai|). (b) with a = 10“^ and (c) with a = 0.1 

are computed for MEG parameter A(fec) = 0.02. Common 
are the exchange parameter D = 0.02, sound velocity ct = 1, 
and sound attenuation rate fp^ = 10“®. Note the change of 
scale between (b) and (c). The insets provide an expanded 
view of the laser spot. 


where Ai(^) = (I/i/tt) sin(fc^) exp(—P/4) and 

A 2 (C) = exp [— -I- l)/(2ct)^] sin(^/ct) and W and 
Og have been by rendered dimensionless as explaind in 
the Supplemental Material [20]. The (purely magnetic) 
first term on the right-hand side represent the exponen¬ 
tial decay of the initially excited wave packet, while the 
second term is a propagating MP mode. The latter de¬ 
cays with the phonon damping rate, hence may have a 
very long mean free path for materials with high acous¬ 
tic quality like YIG (assuming that doping affects the 
magnetization without increasing sound attenuation). 

When the laser spot size is large relative to the MP 
wave length, i.e., Ct <SC 1, according to Eq. (3) the 
ratio between MP amplitude and IFE strength scales 
as Cj A^(fcc), i-e., increases with sound velocity and de¬ 
creases with spot size. In the other limit, ct ^ 1, the am¬ 
plitude of the long-lived signal is inversely proportional 
to Ct, therefore decreases with increasing Ct- We therefore 
estimate this ratio to be maximal e~^^‘^^/TrA'^(kc)/{2a) 
when the laser spot size matches the MP wave length. 
The peak amplitude of MPs in Fig. 2(c) is around 
3x 10“^, in good agreement with e~^^‘^^ynA'^{kc)/{2a) « 
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FIG. 3: (Golor online) Out-of-plane magnetization dynamics 
rriz with Gilbert damping a = 10~'* (dotted curve) and 0.1 
(solid curve) induced by the spot heating by a laser pulse. 
The blue dashed curve shows the displacement profile (Rz)- 
Other parameters are those in Fig. 2. 


2.7 X 10-3. 

Thermal actuation is caused by the shear force 
generated by the laser heating profile FR^{x,t) = 
(C?yCtr/inTi/C'„)0(t)9^ exp(—x^/W^), since the pressure 
wave is decoupled from the spin wave at 0 = 0. The 
asymptotic expression for a ^ A(fcc) becomes 

TO^(a;, t) = mi 

( Ct{A 3 (i) - (l/2)e-‘/"- 

x[A3(ctt-I-i) - A3(ctt - i)]} ’ ‘ 

I ct(l - e““‘cost)A3(i)-^^(-^(^r/d) _ 

[ [A4{cti + x) — A4{cti — x)] ’ * 

( 4 ) 

where A3(^) = ^exp(—^3) A4(^) = exp[—(^3 q, 2 _j_ 

l)/(2ct)3] cos(^/ct). For YIG, the parameter mi = 
'yb(ririinTi/{2MQpctCy) ~ lO^C^Tm-^K”^. Compared 
to Eq. (3), the heat-induced magnetization has (z) odd 
parity in real space, i.e., mz,{—x,t) = —m,z{x,t), {ii) a 
long-lived localized signal near the excitation spot, and 
(Hi) maximum amplitude of propagation at ~ 1, cf. 
Fig. 3. We also plot the amplitude of the thermally gen¬ 
erated phonon wave front that is trailed by the magneti¬ 
zation. 

Two Dimensions — In the following, magnetization is 
oriented along x by an external in-plane magnetic field 
PqH = 50 mT corresponding to 'ypnH/{2 'k) cx 1.4 GHz. 
The spot size W = 1 pm and the dimensionless veloc¬ 
ities are Ct « 0.43 and c; « 0.82. Fig. 4 summarizes 
our main results for the IFE and heat induced dynam¬ 
ics in terms of the out-of-plane magnetization compo¬ 
nent mz ■ We plot a snapshot at t = 5 ns in the x-y 
(film) plane from the calculation with low (a = IQ-"^) 
and enhanced (a = 0.1) magnetic damping in (a) and 
(b), respectively. Fig. 4(a, left) displays IFE actuated 
outgoing rays that broaden with distance from the exci¬ 
tation spot. This feature is insensitive to MFC strength 
and can be understood by the angular dependent group 
velocities of magnetostatic spin wave dispersion around 
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FIG. 4: (Color online) Two dimensional profile of out-of-plane 
magnetization, at t = 5 ns due to (left) IFE field 
(middle) pressure stress Fh, and (right) shear stress ■ We 
normalize the result by mo for (left), m 2 = 'ybrjVFriTiifiX -\- 
2jj.)/{2Mop^CvCti^o) for (middle) and mi for (right). The 
Gilbert damping coefficient a = and 0.1 for (a) and (b), 
respectively. The dashed (olive) curve in (a, left) illustrates 
the angular-dependent spin wave group velocity. 

the average modulus of the wave vectors fco- As discussed 
in the Supplemental Material [20], the group velocity 

Vg ~ §[jfj,oMo/{2ko)] sin29^H/{H + Mq sin^ 0) gener¬ 
ates an expansion of the initial wave packet as shown by 
the dashed (olive) curve, while the the star-like interfer¬ 
ence fringes are governed by the phase velocity. At larger 
magnetic damping, cf. Fig. 4(b, left), the star-like fea¬ 
tures in the x direction are suppressed in favor of MP 
propagation with transverse sound velocity ct- Dotted 
feature around cit are caused by interference of the lon¬ 
gitudinal MP and the damped residue of the initial mag¬ 
netization wave packet with 0 ~ 0, which has relative 
longer lifetime. Note the mirror symmetry with respect 
to the y axis, rriz (x, y) = rriz (—x, y). 

The quadrupolar features in Figs. 4(middle) with 
nodes along the x and y axes and sin 20 symmetry are 
induced by the pressure Fr, caused by a heat pulse. The 
radii of the circular wave fronts correspond to the longitu¬ 
dinal sound velocity. Figs. 4(right) illustrate that a shear 
stress induces MPs that spread with transverse sound ve¬ 
locity Ct,which are thereby clearly distinguishable from 
the pressure induced signals: generates dipolar sym¬ 

metric features with nodes at the y axis, which follows 
from the cos0 symmetry of the MEG coupling. Clearly, 
both heat-induced signals are antisymmetric with respect 
to reflection at the y axis, mz{x, y) = —mz{—x, y), which 
allows discrimination from the IFE response. Moreover, 
we identify a non-propagating signal in the vicinity of 
the excitation spot [see the center of Figs. 4(b, middle) 
and (b, right)], a “smoking gun” for thermally excited 


dynamics. 

It is not easy to predict the absolute and relative mag¬ 
nitude of the two mechanisms for a given light intensity 
from first principles due to uncertainties in the strongly 
non-equilibrium processes after an intense fs light pulse. 
Microscopic theories address the ultra-fast physics of an¬ 
gular momentum and energy transfer from the light to 
the magnetic order [37, 38] and the lattice [34] and should 
ideally be employed to fix the initial conditions for our 
calculations. But also the long-time response depends 
on several temperature and frequency dependent mate¬ 
rials parameters that govern the IFE, light absorption 
etc. Satoh et al. [17] find a Faraday rotation of the 
probe pulse of the order of milli-rad for 110 /im thick 
bismuth-doped iron garnet, which corresponds to a light- 
induced torque of toq ^ 0.004 for a Verdet constant of 
10'^ radm“^T“^ [20]. With thermal expansion coefficient 
y ~ 10“^/K [39], the thermal torque m 2 ^ 4x 10“^(5T/K, 
which can be larger than mg for pulsed laser-induced 
heating [40]. Similar values may be anticipated for mi 
when effective thicknesses C, ~ ym. We should also 
not forget that the fast light-induced demagnetization [1] 
should affect the response directly under the excitation 
spot, but its diffusion should be slower than the ballistic 
response computed here. 

Conclusion and Discussion — We modeled the spa- 
tiotemporal laser-induced magnetization dynamics in 
magnetic thin films, concluding that magnetoelastic cou¬ 
pling is essential for spin angular momentum transport 
because the phonon group velocity is much larger than 
that of the magnons. An experimental study of the sym¬ 
metry of time-resolved magnetization wave fronts radi¬ 
ating from the excitation spot allows discriminating dif¬ 
ferent laser excitation mechanisms, thereby helping to 
answer the long-standing question on the physical ori¬ 
gin of ultra-fast magnetization dynamics, i.e. whether 
it is caused by coherent light-induced magnetic fields or 
sudden heating of the lattice. Moreover, we clarified the 
optimal size of the excitation laser spot to be the MP 
wavelength; for YIG at an applied field of 50 mT it is 
~ 1 /im. 

The essential role of the MEG coupling might have 
larger ramifications. For example, a number of recent ex¬ 
periments on the spin Seebeck effect on YIG came to the 
conclusion that the thermal spin pumping is not caused 
by terhahertz magnons at energies around ksT, but by 
spin waves in a low energy band close to the gap [41-43]. 
Spin information was found to propagate in YIG diffusely 
over large distances [44, 45]. From the present results we 
venture that strongly coupled magnon-polarons could be 
the carriers of thermal spin currents. 

This work is supported by the DFG Priority Program 
1538 SpinGat, the FOM foundation, the E-IMR and IGG- 
IMR, and the JSPS (Grant Nos. 25247056, 25220910, 
26103006). We acknowledge discussions with Benedetta 
Flebus and Rembert Duine. 
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Note added in proof: Ogawa et al. report generation 
of laser-generated magnon polarons that drive magnetic 
bubble domains [46]. 
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This supplement includes the linearized equations of motion in Fourier space with choice of parameters, as well 
as analytical expression for the IFE-induced dynamical response to line and point sources and estimation on the 
magnitude of IFF. 


1. LINEARIZED EQUATIONS OF MOTION IN FOURIER SPACE 

After introducing the forces and torques F acting on the displacement vector = {my,mz, Ri, Rt, Rz)^, the 
linearized equations of motion in Fourier space read [1] 

— 2iujT~^ — (^k^)pRt = ibkcos29my — (SI) 

— 2iujT~^ — Cfk^)pRz = ibk cos Orriz — (S2) 

— 2iujT~^ — c^k^)pRi = ibksu\26my — Fr, , (S3) 

iurriy + ( 7 -ff/j + iauj)mz = i^i cos 9Rz + Fm„, (S4) 

iurriz — {uH'i. + iaoj)my = —iAi(cos 20i?t + s\n29Ri) + (S5) 

where Ai = jbk/Mo, H'f. = Hk RoMq sin^ 9 with Hk = p^H + 7“^(2A/Mo)A:^, and 9 is the angle between magnetic 
field and in-plane wave vector k. a and Tp are the Gilbert damping constant and phonon relaxation time [2]. 

We adopt the material parameters of yttrium iron garnet (YIG). Specifically, the dipolar field and gyromagnetic 
ratio are pqMq = 175 mT and 7/27r = 2.8 x 10^°Hz/T, respectively [3]. The MEG constant b = 6.96 x 10^ J/m^ [4], 
the exchange parameter 2A/Mq = 9.2 x 10“®m^/s [5], and the mass density p = 5170kg/m^ [6]. The transverse and 
longitudinal sound velocities are Cj = 3.8 x 10^ m/s and c; = 7.2 x 10^ m/s [4]. The phonon relaxation rate in YIG 
is of the order of 1 ps~^ [7]. In high quality films, the Gilbert damping coefficient is typically around 10“^ [3]. For 
comparison, we also show results for a = 0.1, which is close to the upper limit for Bi-doped YIG [8]. 


2. ANALYTICAL EXPRESSION FOR IFE-INDUCED ONE-DIMENSIONAL DYNAMICS 


The time evolution of the out-of-plane magnetization Eq. (2) of the main text can be simplified to 

Tno V ' f ,7 /A ^ ( Wq: Of) 


,{x,i) =/ dke ^ /^Im 
, /'TT / I ^ 




(jJzf — i0± 


cos kx, 


where uj± = 


iip + rip zb y (rip — rip)^ + sa^ 


/2. Here Qq = (1 + *«)(! + Dk^), Utp = Ctk -Uifp k = kW, x = xjW, 


fp = UloTp, A = A/Qo ct = Ct/(imo)) i = Ulgt, and D = 2A/(MoQofL^) are all dimensionless. 


3. ANALYSIS OF STAR-LIKE FEATURE DUE TO IFE FROM POINT SOURCES 

The IFE-amplitude induced by an axially symmetric light pulse as plotted in Fig. 4(a,left) displays suprisingly 
complex interference patterns. These features can be explained in terms of the angular-dependent group velocity at 
constant ko = (|k|) [yellow circle in the (volume) spin wave dispersion shown in Fig. Sl(a)]. The thin (black) arrows 

indicate the group velocities Vg = VkWkUo.e — HlToMo/{‘2ko)] sm{29)^F[/{H Mq sin^ 9), where we used the 
simplified dispersion Wk — UPo\jH{H Mq sin^ 9) that is accurate for not too small ko > I/IT. In this approximation 
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FIG. SI: (Color online) (a) Approximate (bulk) spin wave spectrum cuk = 7Mo \/H{H + Mo sin^ 9) in k-space. The thin (black) 
arrows show angular dependence of the group velocities Vg = Vk^Uk at the average modulus fco = (Ikj). Those for three typical 
wave vectors lia,b,c represented by thick (purple) arrows are labeled by Va,b,c- (b) Replotting of Fig.4(a,left), IFE-induced 
magnetization profile in real space at 5 ns, with three thin (black) [thick (purple)] arrows representing tVa,b,c [ka,6,c]. 


Vg • k = 0, i.e., the group velocity is purely transverse. The thick (purple) arrows Fig. Sl(a) represent three typical 
wave vectors ka^b,c, whose corresponding group velocities are labeled by Va,b,c respectively. In Fig. Sl(b), we replot 
the IFE-induced magnetization profile in real space with dashed (olive) curve representing t |vg| and the vectors ka,b,c 
and t'Va^b,c as thick (purple) and thin (black) arrows for t = 5 ns. We can understand this pattern in terms of the time 
dependence of a linear superposition of wave packets arranged around the origin: Each wave packet centered around a 
certain ko is modulated in the direction ko with wavelength l/feg. The whole package moves rigidly in the transverse 
direction with group velocity Vg that is finite by the dipolar interaction. The initial wave packet at the origin thereby 
blows up like |tvg|. The interference pattern initially normal to ko is in real space shifted and combines with the 
patterns generated by other wave packets on the ring to form complex features that in the snapshot plotted in the 
figures look like a star. Moreover, the lifetime of spin waves is determined by I/(awk), leading to Tk„ > Tk,, > Tk„. 
This explains that even for large dampings a radial pattern is still visible near the y-axis while that along the x-axis 
has disappeared [see Fig. 4(b,left)]. 

Analytically, the out-of-plane magnetization at r = {r cos (p,r sin ip) with respect to the center of initial spot can be 
expressed as 

m,ir,ip,t) (X sin" e) cos (^kk cos{0 - p) - . (S6) 


Far from the excitation spot, i.e., r I, the fc-integral is dominated by the spin waves satisfying cos{9 — p) ^ 0 oi 
0 ~ (p ± 7r/2. Taylor-expanding in the phase factor of the cosine leads to 


(r, p, t) - dkke-'^"/^ J dO' cos (k{f - t\vg\)e' - + Mq cos 2 t^)) 


2 cos (tjpoVHiH + Mo cos^ p)'j 


dk 


sin (jt{f — t|ug|)50^ 


-k^/4 


f - t\Vg 


(S7) 


where SO is a small cut-off angle and jugj = (y/roMo/A)] sm2p\y/H/{H + Mq cos^ p). Eq. (S7) catches most features 
observed in Figs. 4(left) as summarized above: (z) cos {t'^pLo'jH{H + Mq cos^ v’)^ is an angle-dependent oscillations, 


i.e., reproduces the star-like feature, {ii) exp[—atjpo\/H{H + MqCOS^ p)] explains the angular dependent attenu¬ 
ation, revealing that rays with p 7^12 decay slower than those for ~ 0. At large damping a radial pattern is 
visible near the y but not in the x-direction [see Fig. 4(b,left)]. {in) The radial integral is a function of r — tjugj only, 
therefore predicts that the wave packet expands with the group velocity jugj. 
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4. MAGNITUDE OF IFE 

Here we estimate the magnitude of the IFE based on the experiments by Satoh et al. [9], who excited spin 
waves in a 110/rm thick Gd 4 / 3 Yb 2 / 3 BiFe 50 i 2 film by laser light with pump fluence of 200 mJ(cm)“^ and pulse 
width of 120 fs at wavelength 1400 nm, and measured Faraday rotations of a linearily polarized probe pulse (at wave 
length 792 nm) of the order of 'dp ~ 1 mrad. Satoh et al. do not measure the Verdet constant. In Bi-YIG (with 
unspecified doping concentration) the Verdet constants are V ~ 5800, 1050, and 380 radm“^T“^ for wavelengths 
543nm, 633nm, and 780nm, respectively [10]. The Verdet constant in (BiYbTb) 3 Fe 50 i 2 at 780 nm was found to be 
V ~ 8.5 X 10^ radm“^T“^ [10]. Here we adopt V = 10"* radm“^T“^ for both pump and probe fields. 

The IFE and FE may be formulated in terms of a phenomenological free energy [11, 12] 

F = Ppnnip{tioM^)Iin, (S8) 

where /3pump is a magneto-optical susceptibility at the pump laser frequency and Lm is the intensity of (right-hand) 
circularly polarized light. The effective field idipp = —{1/^o)9m^F = —Ppumpl'm- In the same theoretical framework, 
the Faraday rotation induced in the external field-free limit hy B = fJ,oMz reads 

= ^^^probe(MoM2)L = Vprohe{f^oMz)L, (S9) 

where Cq, Aq, and L are the velocity of light in vacuum, the wavelength of the probe pulse, and the film thickness, 
respectively. This relation allow us to express (3 in terms of V. 

We can now estimate (assuming Vpump ~ Vprobe) Imo^^ife] = Mol^Ao/in/(27rco) ~ 0.2 T for the experimental 
pump intensity, idipp when applied as a pulse with width r/ generates a magnetization My = yTifioHip^MQ and 
mo = My/Mo ~ 0.004. Disregarding damping, the precession in the external field Hx preserves this amplitude such 
that Mz « My. This leads to r?F ~ V{'yTifilHipEMo)L ~ 0.5 mrad, in reasonable agreement with the observations ]9]. 
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